Experimental magic state distillation for fault-tolerant 
quantum computing 

Alexandre M. Souza 1 , Jingfu Zhang 1 , Colm A. Ryan 1 & Raymond Lafiamme 1 ' 2 

1 Institute for Quantum Computing and Department of Physics and Astronomy, University of Wa- 
terloo, Waterloo, Canada, ONN2L 3G1. 

2 Perimeter Institute for Theoretical Physics, Waterloo, Canada, ON N2J 2W9. 

Any physical quantum device for quantum information processing is subject to errors in 
implementation. In order to be reliable and efficient, quantum computers will need error 
correcting or error avoiding methods. Fault-tolerance achieved through quantum error cor- 
rection will be an integral part of quantum computers. Of the many methods that have 
been discovered to implement it, a highly successful approach has been to use transversal 
gates and specific initial states. A critical element for its implementation is the availability 
of high-fidelity initial states such as |0) and the Magic State. Here we report an experiment, 
performed in a nuclear magnetic resonance (NMR) quantum processor, showing sufficient 
quantum control to improve the fidelity of imperfect initial magic states by distilling five of 
them into one with higher fidelity. 

Quantum information processing (QIP) ^promises a dramatic computational speed-up over 
classical computers for certain problems. In implementation, the physical quantum devices for QIP 
are subject to errors due to the effects of unwanted interactions with the environment or quantum 
control imperfections. In order to be reliable and efficient, quantum computers will need error 
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correcting or error avoiding methods. One method to achieve fault-tolerant quantum computation is 
to encode the state of a single quantum bit (qubit) into blocks of several qubits that are more robust 
to errors. Based on this idea, quantum error correction codes, the theory of fault-tolerant quantum 
computation and the accuracy threshold theorem have been developed 5 1 . A key element for fault- 
tolerant quantum computation is to avoid bad error propagation. One straightforward protocol is to 
use transversal gates where an error occurring on the A;th qubit in one block can only propagate to 
the kih. qubit in the other blocks. A highly successful approach to achieve fault-tolerant universal 
quantum computation is based on quantum error correcting codes with gates from the Clifford 
group that can be applied trans versally^. Unfortunately they are not universal l±21JJJ and they must 
be supplemented with the preparation of not only the |0) state but also another type of state such 
as a. Magic State^^. Thus, a critical element for fault-tolerance is the availability of high-fidelity 
magic states. Consequently, in the pursuit of experimental fault-tolerant quantum computation, it 
is important to determine whether we have sufficient experimental control to prepare these magic 
states. In general these will be prepared with some imprecision. The states can be improved 
by distilling many magic states to produce a fewer number of them which have higher fidelity. 
Here we report an experiment, performed in a seven-qubit nuclear magnetic resonance (NMR) 
quantum processor, showing sufficient quantum control to implement a distillation protocol based 
on the five-bit quantum error correcting code ^^21 which uses only Clifford gates. The fidelity of 
imperfect initial magic states is improved by distilling five of them into one with higher fidelity. 
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Results 



Theoretical protocol. The Clifford group is defined as the group of operators that maps the Pauli 
group onto itself under conjugation. The Pauli group is defined as 1 {±11, ±il, ±<r x , ±icr x , ±cr y , 
± io-y, ±a z , ±ia z } where a x , a y , a z , and 1 denote the Pauli matrices and identity operator, respec- 
tively. The Clifford group on n qubits is a finite subgroup of the unitary group U(2 n ) and can be 
generated by the Hadamard gate H, the phase-shift gate S p h, and the controlled-not gate CNOT 
represented as 
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An arbitrary one-qubit state can be represented in the Bloch sphere as 

p = (l + p x a x +p y (Ty +p z a z ) /2 (2) 

where p x , p y and p z are the three polarization components of the state. The magic states ^ are 
defined as the 8 states with p x = ±1/ \/3, p y = ±1/ a/3, p z = ±1/ v3 (T type) and the 12 states 
with p x = 0, p y = ±1/V2, Pz = ±1/V2; Py = 0, Pz = ±1/V2, p x = ±1/V2; p z = 0, 
p x = ±1/V2, p y = ±l/y/2 (H type). These states are called "magic" because of their ability, 
with Clifford gates, to enable universal quantum computation and the ability to be purified, when 
it has been prepared imperfectly, using only Clifford group operations ™ In our current work we 
distill an imperfect magic state into a T-type magic state represented as 

p M =[l+(a x + a y + a z )/V3]/2. (3) 
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To quantify how near a state p is to the magic state, we define the m-polarization (polarization in 
the direction of the magic state) 

p = 2 Tr[p M p] - 1 = -j= {p x +P y + p z ) ■ (4) 

The distillation algorithm requires five copies of a faulty magic state p in = p m as the input 
state. In the original proposal 12 , the measurement of four stabilizers Si (i = 1, 4) is applied to 
Pin, where Si = a x ® a z ® a z ® a x ®% S 2 = 1 <8> cr x <g> o z <g> o z ® a x , S 3 = <j x <g> 11 <g> a x ® a z <g> a z 
and S4 = o z ® a x ® 11 ® cr x <g> o z . If the outcome of any of these observables is —1, the state 
is discarded and the distillation fails. If the results of all the measurement are +1, corresponding 
to the trivial syndrome, one applies the decoding transformation for the five-qubit error correcting 
code 17 to the measured state and obtains the output state p dis <g> 1 0000) (0000 1 where p dis has the 
output m-polarization p out . If the input m-polarization p in > p = w 0.655, distillation is 

possible and p out > Pi n and produces a state nearer to the magic one. In an iterative manner, it is 
possible to obtain the output m-polarization approaching 1. 

As NMR QIP is implemented in an ensemble of spin systems, only the output of expectation 
values of ensemble measurements ^ are available. Consequently the above projective measure- 
ment of the stabilizers cannot be implemented in our experiment. However, as the decoding opera- 
tion is just a basis transformation from one stabilizer subspace to another, it is possible to evaluate 
the result of the distillation after decoding. Therefore, we directly apply the decoding operation 
to the input state p in , and the output state becomes a statistical mixture of 16 possible outcomes 
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represented as 

15 

Pout = J2 d iPi ® ( 5 ) 

where ^ is the probability of each outcome, and \i) = |0000), |0001), ... , 1 1 1 11), for i = 0, 
1, 2, 15, noting p = pdis- Now measuring |0) on all four qubits in \i) indicates a successful 
purification. We can obtain 9i and p t using partial quantum state tomography^. 

Experimental results. The data were taken with a Bruker 700 MHz spectrometer. We 
choose 13 C- labelled trans-crotonic acid dissolved in d6-acetone as a seven-qubit register. The 
structure of the molecule and the parameters of qubits are shown in Tabled! We prepare a labelled 
pseudo-pure state p s = 00cr z 0000 using the method in Ref. where = |0)(0| and the order of 
qubits is M, Hi, H 2 , Ci, C2, C3, C4. One should note that we are using the deviation density matrix 
formalism. 

We prepare an initial imperfect magic state with three equal polarization components by 
depolarizing the state = (t + a z )/2. First we apply an/2 pulse to rotate the state to (1 + a x ) /2 
and then another ir/2 pulse along direction [cos a, sin a, 0] is applied. We use phase cycling to 
average the x- and y- components of the state to zero, and therefore the polarization of the spin 
initially in the state is reduced. The depolarized state is represented as 

(l-<7,sina)/2. (6) 

Finally we apply a rotation with angle arccos(l/A/3) about the direction [— 1/a/2, 1/a/2, 0] to 
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obtain an imperfect magic state 

p=[l + p(a x + o y + a z )/V3] /2 (7) 

where p = — sin a. The evolution of a z in the preparing p is shown in Fig. [TJ By doing the above 
operation for qubits M, Ci, C 2 , C 3 , C 4 , respectively, we obtain five copies of the imperfect magic 
states pi n = p® 5 . Exploiting partial state tomography, we measure p for each qubit and use the 
average as the input m-polarization p in for p in . 

The circuit for the distillation operation is shown in Fig. |2] Ci carries the distilled state 
after the distillation. With partial state tomography, we can determine 6i and pi in equation ©, 
where p = Pdis, from which the output m-polarization p out is obtained. The experimental results 
for magic state distillation for various p in are shown in Fig. [3] where figures (a) and (b) show the 
measured p out and Q Q , respectively. The straight line in figure (a) represents the function p out = Pi n . 
The data points above the line show the states that have been distilled experimentally. 

The implementation time of the distillation procedure is about 0. 1 s, a non-negligible amount 
of time (10%) compared to coherence time (T 2 in Table [T]). Hence the decay of the signals due 
to the limitation of coherence time is an important source of errors. We extract p out by measuring 
the ratio of 9 p out and 9 , where these two factors are obtained by various single coherent terms in 
a series of experiments (see details in Methods). We have assumed that the terms have the same 
amount of decoherence. The results of simulations with dephasing rates T 2 * and T 2 are shown in 
Fig. [3]as the blue squares and red triangles. The simulation results show that the decoherence rates 
are long enough to allow the distillation and suggest that the deviation of #0 from the theoretical 



expectation can be mainly attributed to relaxation effects. Additionally, imperfection in the shaped 
pulses and inhomogeneities of magnetic fields also contribute to errors. 

Discussion 

We modify the original distillation protocol by avoiding the projective measurement, which is 
not possible to implement in the NMR QIP's. We exploit partial-state tomography to obtain each 
output in the mixture of the outcomes after the distillation, and only in a post-processing step do we 
choose the one we need. Although we could access the |0000) (0000 1 subspace using a procedure 
similar to the pseudo-pure state preparation, the method would take a substantially longer amount 
of time and would be more error-prone. In this work we aim for a quantitative result, i.e. increasing 
the magic state purity. We need to minimize the readout manipulations to avoid control-error 
induced distortions of the inferred final state and associated purity. Hence, we limited ourselves to 
simple high-fidelity readout procedures. 

In summary, we have implemented a protocol for distilling magic states based on the five 
qubit quantum error correction code. We exploit five qubits by controlling a seven-qubit NMR 
quantum information processor. The experiment shows that we have obtained enough control to 
purify faulty magic states through distillation. 

Methods 

Overview of the experiment. To implement the experiment, we exploit standard Isech and 
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Hermite-shaped pulses as well as numerically optimized GRAPE pulses to implement single- 
spin operations. The GRAPE pulses are optimized to be robust to radio frequency (r.f.) inho- 
mogeneities and chemical shift variations. All pulses are combined in a custom-built software 
compiler 22 . The compiler loads the information about the internal Hamiltonian and the desired 
unitary transformation from simple predefined building blocks. The blocks are then systemati- 
cally put together to form a pulse sequence ensuring that the errors in the building blocks do not 
propagate as the sequence progresses. 

R.f. selection. The effect of pulse imperfections due to r.f. inhomogeneities is reduced 
by selecting signal based on r.f. power The signal selection is achieved by spatially selecting 
molecules from a small region in the sample. The method is similar to imaging methods^ and has 
been used in previous works 20 . Here we substitute the original pulse sequence proposed in RefP^ 
by a single GRAPE pulse to optimize the performance. Besides reducing r.f inhomogeneities, the 
spatial selection of spins can also reduce the static field inhomogeneities and therefore reduces the 
loss of signal during the experiment. We have found that the effective relaxation time (T 2 *) of spins 
after the r.f. selection increases significantly, e.g., up to a factor 2 for some spins. 

Partial state tomography. We use the spectra obtained from the labelled pseudo-pure state 
p s = 00(j 2 0000 shown in Fig. @]as a phase reference and to normalize the signals in Ci and C2 
spectra, for measuring the initial and output m-polarization. To obtain the reduced density matrix 
of Ci through the partial state tomography, we expand equation © as a sum of product operators 
and represent pi as 

Pi = ^ i 1 + Pi,x&x + Pi,yOy + Pi,z^z) ■ (8) 
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In the expansion there are 128 terms that are required to determine by the experiment. 

The coefficients of such expansion can be directly related to the measurable spectral ampli- 
tudes 19 . On the other hand, such coefficients can also be related to the relevant parameters of ©, 
i.e, Pi ;X , p ijV , p i:Z and 6>j for i — 0, 1, 2, 15. The relation between these parameters and the NMR 
observables can be expressed by the set of linear equations 

C = A x R. (9) 

The nth element, C (n) , of the column vector C is the coefficient related to the operator er™ 4 looser™ 1 a™ 2 a™ 
with the order of qubits M, Hi, H 2 , Ci, C2, C3, C4, where s can be one element of the Pauli group 
{a x , a y , a z , 1} and the vector n = ("hi, n 2 , ^3, ^4) is the four digit binary representation of the 
integer n - 1 . For s = a x , a y , a z , and J, R(n) = 9 n p n ,x, 9 n Pn, y , 0nPn,z and 9 n , respectively. The 
elements of the matrix A are given by 

A(k,m) =Uj =1 (-l) hmi . (10) 

Providing that we have all necessary coefficients measured, we can reconstruct the distilled 
states using the following approach. First, we fit the NMR spectral lines to the yield complex 
amplitudes for measuring all necessary coefficients Fig. \5\ illustrates the spectra of Ci after 
the completion of distillation for p in = 0.95, where the experimentally measured, fitted, and ideal 
spectra are shown as the red, black, and blue curves, respectively. Then the state © is reconstructed 
by solving the set of equations ©. Our calculation shows that four readout operations are sufficient 
to determine all coefficients: first, read out on Cj; second, read out on Ci after the application of a 
7r/2 pulse; third, read out on C2 after the application of a 7r/2 pulse; and forth, read out on C2 after 
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a polarization transfer from Hx to C 2 . The last two readout operations are sufficient to measure 
all 9i, and the first two are used to measure 9^. The errors for the coefficients, as well errors for 
Pi,x, Pi, y , Pi,z and Oi , are estimated from the uncertainty of the fitting parameters. The measured 
initial and output m-polarization, as well as Q { and p h are listed as Supplementary Tables S1-S10 
in the Supplementary Information 25 . The comparison of the various measured p with the theory 
is shown as equations (1-7) in the Supplementary Methods in the Supplementary Information. 
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Table 1 : Characteristics of the molecule of trans-crotonic acid. Molecular structure is shown 
as Fig. [U The chemical shifts and J-coupling constants (in Hz) are on and below the 
diagonal in the table, respectively. The transversal relaxation times T 2 measured by a 
Hahn echo and T 2 * calculated by measuring the width of the peaks through fitting the 
spectra are listed at the bottom. The chemical shifts are given with respect to reference 
frequencies of 700.13 MHz (protons) and 176.05 MHz (carbons). The molecule contains 
nine weakly coupled spin half nuclei but consists of a seven qubit system since the methyl 
group can be treated as a single qubit using a gradient-based subspace selection^. 
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Figure 1: The evolution of o z in the preparation of a faulty magic state in Bloch sphere, 
where arrows represent the states of the qubit. a, A 7r/2 rotation along y-axis transforms a z 
(blue) to a x (green), b, Another n/2 rotation along direction [cos a, sin a, 0] transforms a x (green) 
to o x cos 2 a + a y cos a sin a — a z sin a (black). In phase cycling we apply the second n/2 rotation 
by changing a to it + a to transform a x to a x cos 2 a + a y cos a sin a + a z sin a. After averaging 
the x- and y- components to zero, the polarization is reduced to — a z sin a, shown as the yellow 
arrow in figure c, noting that a G [n,3n/2]. d, A final rotation with angle arccos(l/v / 3) along 
[-1/V2, l/y/2, 0] transforms — a z sin a (yellow) to — sin a(a x + a y + a z ) / y/3 (red). 
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Figure 2: Gate sequence for Magic State distillation. The sequence is constructed from the five 
qubit error correction code ^ where ±X = exp(=Fi7r<T :r /4), ±Y = exp^inay/i), and Z a = 
exp(—iaa z /2). Qubits labelled as M, Ci, C 2 , C 3 and C 4 are used to encode the five copies of the 
initial state. Due to the nature of the algorithm, the carbon Ci contains the distilled magic state 
only when M, C 2 , C 3 and C 4 are in the |0000) state. It is important to emphasize that all gates are 
Clifford gates. The refocussing pulses (which also decouple H\ and H 2 ) are not shown. 
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Figure 3: Experimental results after the completion of magic state distillation. Output m- 
polarization of the faulty magic state (a) and the probability 9 (b) of finding this state in the mix- 
ture of outcomes [see equation ©] as a function of the input m-polarization of the initial faulty 
magic state. The experimental data are represented by the filled circles and the error bars are esti- 
mated from the uncertainty of the fitting parameters. The line in figure (a) represents the function 
Pout = Pin- The experimental points above the line show the states that have been distilled, while 
the points below the line show the states that cannot be distilled in the protocol. The theoretical 
prediction is represented by the black solid curves. The blue squares and red triangles, connected 
by dashed lines for visual convenience, are the simulation results where the dephasing rates are 
chosen as T 2 * and T 2 (see Tabled]), respectively. The effective T 2 during the experiment should 
be similar to the Hahn echo T 2 . The deviation can be attributed to other error sources (see text). 

The dephasing times of Hi and H 2 actually do not influence the results because Hi and H 2 can be 

17 

effectively assumed in and a z during the whole experiment, respectively. 
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Figure 4: NMR spectra. a,b, Spectra of Ci and C 2 obtained by 7r/2 readout pulses when the 
system lies in the labelled pseudo-pure state p s = OOo^OOOO. The vertical axes have arbitrary but 
the same units. 
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Figure 5: Spectra of Ci after the completion of distillation for p in = 0.95. The spectra are 
divided in four different parts shown as figures (a-d) for better visualization. The vertical axes 
have arbitrary but the same units. The experimentally measured, fitting, and ideal spectra are 
shown as the red, black, and blue curves, respectively. 
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Figure 6: Molecular structure of trans-crotonic acid. 
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